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ПРИМЕНЕНИЕ УРАВНЕНИЙ НЕЙТРАЛЬНОГО РАВНОВЕСИЯ 
К ЗАДАЧЕ О ЦИЛИНДРЕ, СЖАТОМ БОКОВЫМ ДАВЛЕНИЕМ 


Г.И. ВОЛОКИТИН, К.В. КОЗЛОВ, В.Б. ФЕДОСЕЕВ 
(Донской государственный технический университет) 


Рассматривается задача потери устойчивости полого кругового цилиндра. Используется теория малых де- 
формаций, наложенных на конечную. Докритическое начальное напряженное состояние неоднородное. За- 
кон состояния определен соотношением Блейтц-Ко. 

Ключевые слова: конечная деформация, напряжения, круговой цилиндр, тензор, выпучивание, устойчи- 
ВОСТЬ. 


Введение. Согласно феноменологическому подходу, принятому в механике сплошных сред, ба- 
зовые соотношения для нелинейно-упругих тел включают конкретное представление удельной 
потенциальной энергии деформации. Предлагаемое выражение должно допускать описание ос- 
новных нелинейных эффектов и приводить к доступному математическому рассмотрению простых 
деформаций. На свойства резиноподобных материалов сильно влияют химический состав и про- 
цессы изготовления. Для некоторых сортов резины подходит модель несжимаемого упругого ма- 
териала, конкретизированная потенциалами Муни, Трелоара, Г.М. Бартенева и Т.Н. Хазановича, 
К.Ф. Черных и И.М. Шубиной [1]. Однако часто необходим учет сжимаемости материала. Как 
практический, так и теоретический интерес в этом случае представляет потенциал Блейтца и Ко. 
Целью данной работы является исследование устойчивости равновесия толстостенного цилиндра 
из этого материала при гидростатическом сжатии на основе трехмерных уравнений нелинейной 
теории упругости. 

Определяющие соотношения и начальная деформация нелинейно упругого цилиндра. 
Закон состояния задается выражением плотности потенциальной энергии деформации, предло- 
женным в работах Блейтца и Ко [2, 3]: 





иВ а ы а 
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Здесь [, [,, [, -— главные инварианты меры деформации Коши-Грина, коэффициенты 


о, В, и - модули упругости. Константа а выражается через постоянные Ламе и коэффициент 


Пуассона: 

А, У 

а = — = . 

2 1-2 
Трехконстантный закон состояния (1) удачно подходит для описания упругих свойств некоторых 
сортов резины, требующих учета сжимаемости. 

Далее будет использован закон состояния для напряжений в форме Фингера [1]: 
2 О О О ОИ 
= р +1 Е Е” 
в а, ог, дог, 
где Т - тензор напряжений Коши; Е - единичный тензор; Р’- мера деформации Фингера. 
Выражение тензора Т , найденное с учетом потенциала (1), принимает вид: 


п [сПв+еЕк-с'Е? |. (2) 
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Обобщенные модули задаются равенствами: 
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Рассматривается осесимметричная деформация кругового полого цилиндра. 


27 


Технические науки 








Материальные координаты — цилиндрические координаты точки недеформированного цилиндра -— 
’, ф, 2. Координаты точки после деформации, возникающей за счет поверхностного нагруже- 
ния равномерно распределенным давлением по наружной боковой поверхностям цилиндра, обо- 
значены соответственно К, Ф, 7. Место точки в актуальной конфигурации определяется вы- 
ражениями: 

К=К(г)е, +4, Ф=ф. (3) 
Здесь е., е‚, &,- базисные векторы цилиндрических координат; г, ф, 2 -радиус, азимуталь- 


ный угол и высота соответственно. Введем обозначения: 
К 
а= В', р=-—. (4) 
Г 
Тогда градиент деформации, определяемой соотношениями (3), принимает вид 
о ото д ой 
УА= | и —ё, и 1 —- К=ае,е, + ее, + ЧЬ . 
и г (©) 7 
0 
где У - набла-оператор в отсчетной конфигурации. Отсюда находим меру Фингера и главные 
инварианты: 


Е=а’ее, +Ь?е,е, +4, П=1(Е)=а’ +? +а?, 1, =а? (6? +4?)+574?, 1, =а’?4?. (5) 
Тензор напряжений является соосным мере Фингера и с помощью физических компонент [1] вы- 


ражается в виде: 


Т=ос,ее, +0,е.е, +6 В, 


где, учитывая формулы (2) и (3), для физических компонент можно получить следующие выра- 
жения: 


с, = М НУ ® а? (ур, +1)-а“), 
ь О НУ“ -+Ь (у1,+1)-Ь*), 
с. = ь НУ 1, +4? (м,+1)-4“). (6) 
В соотношениях (6) инварианты вычислены по формулам (5), а у, = т. 


Если отсутствуют массовые силы, условие равновесия в объеме имеет вид 
АА ВЕ с О. ы 
+ +1 (нее, +06, +0.4& ) =0. 





е, е 5 

ОК К 'дФ 07 
Учитывая зависимости (6) и деривационные формулы [1], можно это условие свести к одному 
уравнению 





4с, 
= +5 (6,-0,)=0. (7) 


Поскольку в этом уравнении напряжения с., с, определены формулами (6), то ясно, что (7) — 
это нелинейное дифференциальное уравнение второго порядка. Приведем вариант уравнения, 





отвечающий упрощенному варианту потенциала (1), - материалу ’Ноулса-Стернберга 
(=... В=0;. 14): 
= 2 3 
Ча Ь А ый (8) 
Я З1ЦЬ ВБ Ь 
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Условия равновесия на цилиндрических поверхностях имеют вид: 


0, ГЕК 
о, = ; (9) 


_Р, ГЫТ, 
где = внутренний радиус трубы; р = наружный радиус; р - интенсивность давления. В ча- 


нН 


1 
стности, для материала Ноулса-Стернберга с, = ,(1- ==] 
а 
Краевое условие на торце удовлетворяется интегрально, в «смысле Сен-Венана»: 
6 (т) 
О= | с.абгаГ, (10) 
вЫ) 
где О - продольная сила, приложенная к торцам цилиндра для равновесия. Окончательно зада- 
ча Ламе для цилиндра определяется краевой задачей (7), (9). Причем параметр осевого смеще- 
ния точек 4 задается в соответствии с условием (10). В упрощенном случае материала Ноулса- 
Стернберга задача Ламе для нелинейно-упругого цилиндра рассмотрена в [4]. 

Для решения нелинейной краевой задачи (7), (9) предложен численный метод, исполь- 
зующий конечные разности [5]. Выбираются равноотстоящие узлы, производные аппроксимиро- 
вались центральными разностями. Узловые значения аргумента определяются формулами: 

$ р ЕЕ 
п = +(/-2)1, 1=2,3,...М+2, в 
Используем также еще две «законтурные» точки и =и,-й, гу. =7, +й. Обозначим х, значения 


искомой функции К(ь) на выбранной сетке значений аргумента. Тогда, в частности, из (7) имеем 
систему № +1 нелинейных алгебраических уравнений: 








р Г В. нь 5 
А АЯ 2 з |! 
2 Эй, Ь ББ Ь 
Хы НХ, х, . 
где а= 5» —, р=—^ . Еще два условия определим из краевых условий: 


и. 


—2/3 —3 —2/3 
Хх 
вю) | пота [5=) | 
2 н 


Таким образом, имеем систему №-+3 алгебраических уравнений с М+3 неизвестными. Решение 
системы получено методом простых итераций [5]. В качестве нулевого приближения итерацион- 
ного метода решения краевой задачи на сетке принимались узловые значения соответствующего 
решения задачи Ламе для цилиндра в линейной теории упругости. Такой выбор обусловлен тем, 
что при малой относительной толщине оболочки и небольших давлениях уравнения (7), (8) и (9) 
переходят в соотношения линейной теории, для которых известны точные решения [6]. 

Следует отметить, что в упрощенном варианте материала Ноулса-Стернберга допускается 
решение дифференциального уравнения (8) в аналитическом виде. Его можно переписать, учи- 











а а-Ба 
тывая тождество — = —— ‚ следующим образом 
а’ г @ 
аа Па а’ а 
+ +3 |=9 
а 6 Ь Ь 
а аа И Ч а-Ь 
Применим подстановку #=—. С помощью тождеств — =#+Ь и = последнее диффе- 
Ь аь ар 4 г 


ренциальное уравнение можно представить в виде: 
Аг Ге 


к ЕТ 
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Отсюда находим 








3 
|4 3 2+1 
ЕС, ый агсфап = Слу(Ю, (11) 
1 т / <] 

(+1448 “М ы 

где С, - произвольная постоянная. 
@ ТАК 
Далее, учитывая равенство — = т получаем другое дифференциальное уравнение 
| 
АК @ 


ОЙ реа. 





Отсюда имеем 








1 
(2 +1+4)4 3 ат 
В=С ———__——ех атсфап =С.6(0. (12) 
д) СК 


Выражения (11) и (12) дают общее решение исходного дифференциального уравнения (8) 
в параметрической форме. Для решения краевой задачи необходимо удовлетворить краевым ус- 
ловиям (9) и найти границы изменения параметра #, и Е (В н>и, Е н>и, ). Эти значения пара- 


метра г определяются из системы двух нелинейных уравнений, вытекающих из краевых условий (9): 
п, _ У(ь) 
Г (1 )' 
ы 
|) р 7, Са) 
| =|1+— | — : 
1. | т. е (+. ) 


Далее, изменяя параметр 1 от 1, до (,, с помощью (11) и (12) можно определить зависи- 





= Ш 





мость А(г). В [4] указаны полученные ранее другими авторами параметрические формулы обще- 
го решения дифференциального уравнения (8). Однако соотношения (11), (12) имеют отличия, 
связанные с выбором параметра. Также в [4] не прояснен вопрос решения краевой задачи. Точ- 
ное решение сходной нелинейной краевой задачи Ламе для сферы из материала Ноулса- 
Стернберга приведено в [7]. 
Уравнения нейтрального равновесия. Вектор добавочного перемещения возьмем в виде, до- 
пускающем несимметричные формы выпучивания: 
и=и(г,ф,2)е, +у(г,ф,2)е, + иг, ф,2)ь . 

Применяя деривационные формулы, для градиента вектора добавочной деформации име- 

ем равенство 


1 ди 1 ду 1 д 1 ду 1[ ды 
Уве и еее ие 
а д" К [610] Я 02 а ди К\ д 


1 ди. 1 д». [0 о 10. 
це, + Бе, + е.1; + ев. (13) 
а 82 а 02 а 0’ К дф 
Смежные формы равновесия цилиндра исследуем с помощью оператора А.И. Лурье 
© [1,6], который для материала Блейтца и Ко определен соотношением (10), полученным в ра- 
боте [7]. Учитывая это соотношение и равенство (13), для тензора ® имеем следующее пред- 


ставление: 
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9 = ИА ве и е„е, + В вв; а ый № ее, + 
д 0Ф де 0Ф В 


02 74 


+ Е. а ем + ОТВ ЗВ е„е + 
д 0Ф 02 д 0Ф ы 


+ ПОВ И ее Ои- И 
ди Оф ди 02 


ди ди \. ди д} . ди ду]. 
+ а 


Коэффициенты при неизвестных и, у, м’ и их частных производных известны после решения 


(14) 


начальной задачи Ламе функции переменной 7. Приведем выражение одного из этих коэффици- 
ентов: 





1- 
а а ((1 +24) 17° + (2а-1) “9 + у а?) + 35? 4?) 
3 


Уравнения нейтрального равновесия в объеме [1] (У® =0) приводят к системе трех ли- 
нейных дифференциальных уравнений второго порядка с переменными коэффициентами: 


09 
р. ты 9,9, +—— |+ ео 0, 
В д ] а д 


==) а 00., —_ 





00, а 
— ве. НО 





+ =0, 15 
6 а @ и 


6) 
00 “о. = 229 


д’ 


0. 





07 К 








= 
ОФ Я 02 

В уравнениях системы (15) компоненты тензора ® выражены в соответствии с представлением 
(14). Нейтральное равновесие на боковой поверхности определяется условием [1,6] 

№ =-р(Узм- Уи”), (16) 
где У - набла-оператор в метрике деформированного состояния; № - вектор внешней нормали 
к деформированной поверхности, нагруженной гидростатическим давлением интенсивности р. 
На торцах принимаем граничные условия вида при 2=0 и 2=[ ( Г- длина цилиндра): 

и=0, ©9_=®©. =0. (17) 
Соотношения (17) означают недопустимость добавочного перемещения в осевом направлении и 
отсутствие трения на торцах [8]. Эти равенства не противоречат характеру нагружения цилиндра 


наружным давлением и сжатия с торцов абсолютно жесткими гладкими плитами. 


Разыскиваем решение краевой задачи (15), (16) 
и=Х „(г)с0зифсо$ 2, 
У=У„(Г)5шифсо5 А, „2, (18) 
"=й„(г)соз ифзшА2, 


пт 


ге пи ^, - параметры волнообразования, в соответствии с (17) определяется равенством 


пт 
А ва - числа п, т -— целые. 


т 
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Подстановка (18) в уравнения (15) и краевые условия (16) приводит к разделению пере- 
менных. В результате получаем однородную систему обыкновенных дифференциальных уравне- 
ний: 


АХ", +| 4+4 +5 (4 —В, )|х. +4 Нави + |4 а, > + 





+| 4+4 В, -т?Е, ) м, Хин ще (-вЬ- В) |. 


я | + (4 В: )|2. 20 


пт 


ру" + Б!-р,+=(Р.+Е) |", -п| р, +В |Х,-| + (р, +8, + 
р: |: р: 
БРВ ОДИ,. п рр. +В,) | Х.- (19) 


а а 
—п^ и [18 + 50. == 0, 


пт 


С", +| +в, Р.А, 6+5, Х-, ба, +5 С, Х„- 
р. а К а 
В +4, 1. = 
у а 


пт 


[ая с, =0 


Из краевых условий (16) после разделения переменных получим при и=и, иг=и: 


пт 


АХ! +4 о +7) + [4 Ч . = 0, 


пт 


ру, -[2.-$)ох +7„)=0, (20) 


пт 


А. [с 1 = 0. 
Ч 


В условиях (20) 4=0, если г=х",, и а=р, если г=и.. 
Однородная система (19) и краевые условия (20) выражают нейтральное равновесие ци- 
линдра. При бифуркационных значениях параметров р и 4 возможны нетривиальные решения 


краевой задачи (19), (20). Они соответствуют возмущенным равновесным формам цилиндра. Кри- 
тические значения давления находим как наименьшие бифуркационные значения р при надле- 


жащем выборе и и /.„. Бифуркационные значения — собственные числа задачи (19), (20). Отме- 
тим, задача на собственные значения является нелинейной, поскольку искомый параметр входит 
нелинейно через функцию В(,). Для решения такой задачи на собственные значения использо- 


ван численный метод [8]. 
Выводы. Численный алгоритм определения бифуркационных нагрузок, включающий решение 
задачи Ламе о начальной деформации цилиндра, вычисление коэффициентов системы (19) и по- 
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лучение собственных значений, реализован на ЭВМ. Расчеты проводились как для упрощенного 
варианта материала — материала Ноулса-Стернберга, так и для полного закона Блейтца и Ко. В 
качестве примера применения полученных уравнений нейтрального равновесия рассматривались 
некоторые случаи потери устойчивости, изучаемые в теории оболочек [9]. Для очень длинных 
цилиндрических оболочек под внешним боковым давлением критическая нагрузка разыскивалась 


в виде соотношения Грасгофа-Бресса: р=2н(1+у) р.=”, где = - относительная толщина обо- 
лочки, р. — безразмерный параметр давления. Обнаружено, что критическому (минимальному) 


давлению отвечают плоские формы выпучивания: и=2, т=О. Представления о влиянии нели- 


нейности на величину верхнего критического давления в зависимости от относительной толщины 
отражены в табл. 1. Данные второй строки согласуются с результатами, приведенными в [8]. Чис- 
ленный анализ также показал, что для не слишком толстостенных цилиндров влияние на крити- 
ческое значение параметра давления точности решения начальной задачи (задачи Ламе) и кон- 


В 


К К И. 
станты В незначительное. Изучался случай оболочек средней длины, где принималось = 0,38. 


Для таких оболочек критическая нагрузка разыскивалась в соответствии с формулой Саусвелла- 
Папковича: р = 2и(1+у) р.=”. Длина оболочки неизменная (4 =1), цилиндр нагружен внешним 


боковым давлением. Критические значения параметра давления и соответствующие значения па- 
раметров волнообразования содержатся в табл. 2, параметр т=1. 









































Таблица 1 
&=2(»-*,)/(". +") 0,05 0,10 0,20 
р. (у=0,272) 0,26 0,27 0,29 
р. (у=0,250) 0,34 0,36 0,38 
Таблица 2 
&=2(н,-—^,)/(",+›,) 0,05 0,10 0,20 
Й 4 Е 3 
р. (У=0,272) 0,38 0,39 0,44 
р. (у=0,250) 0,49 0,52 0,56 
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АРРЫСАТТОМ ОЕ МЕЧТКАЕ ЕОУТЫВВТОМ ЕОЧАТТОМ$ 
ТО РКОВЕЕМ ОМ СУММОЕК М/ТТН ГАТЕКАЕ РКЕЗЗУВЕ СОМРКЕ$З$ЗТОМ 


С.Т. УОБОКИИТМ, К.\У. КО2ТОМ, М.В. РЕРОЗЕУЕМУ 
(Роп Зае Теспткса! Упмег®Ку) 


Тре БисЮта ргоМет о! пе по!ои/ сисшаг сУпаег 15 сопаегед. Тре зта/-5гат {теогу зиреттрозей оп те 
ИпКе 15 изеа. Тре зибспйса! па! 5#тез$ 5Е ие 15 протодепеоиз. ТПе [аи/ оЕ ее [5 аейтеа Бу В/а-Ко. 
Кеуигога5: Ппйе цеюгтайоп, 5{те5е5, сгсшаг су[ПаЕЕ, {епзог, Бискта, 5аИКУ. 
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